Viscosity of a two-dimensional Fermi liquid 
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Shear viscosity of a two-dimensional Fermi liquid is found to be a nonanalytic function of temperature. In 
contrast to the quasiparticle lifetime that is determined by the forward-scattering processes, the main contribu- 
tion to the viscosity arises from the quasiparticle scattering in the Cooper channel. The viscosity is enhanced by 
the logarithmic singularity of the vertex part. This singular behavior can manifest itself in the two-dimensional 
electron transport, and in the momentum relaxation of fermions in atomic traps. 
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Viscosity characterizes momentum relaxation in a fluid. In 
the presence of a weakly inhomogeneous flow u(r), the linear 
relation between the stress tensor of the d-dimensional fluid 
and the flow gradients is conventionally written as [1] 



(1) 



Iljj = -?? diUj + djUi - -*, 3 Vu - CVu 



The coefficients 77 and £ are the first (shear) and the second 
(bulk) viscosities, out of which 77 usually dominates: 77 3> (. 
The linear relation (1), complemented by the momentum con- 
servation law p (pt + uV) Uj = —diHij ~ rjS? 2 Uj, where 
p(r) = mN(r) is the fluid mass density, leads to the Navier- 
Stokes description on the length scales exceeding the molec- 
ular mean free path t. Such a hydrodynamic approach bears 
great predictive power, since a single parameter 77 accounts for 
the details of the molecular interactions on the scale I. The 
dependence of the viscosity on the temperature and density 
directly determines the fluid's relaxation and transport prop- 
erties. Calculating rj(T, N) is in general a difficult task. 

Early developments of the kinetic theory by Chapman and 
Enskog have lead to the understanding of the viscosity of di- 
lute gases. The latter can be estimated as [1, 2] 



77 ~ pvi , 



(2) 



where v is the typical velocity of a molecule in a gas, and the 
mean free path £ is determined by the molecular density and 
the scattering cross-section. Assuming weak energy depen- 
dence of the scattering, one readily obtains the universal tem- 
perature dependence of the gas viscosity 77 oc VT that stems 
from that of v according to the Maxwell distribution. 

Later on, the notion of viscosity was generalized for hydro- 
dynamic modes in quantum fluids [3], In particular, the shear 
viscosity has been studied in the context of 3d Fermi liquids 
(FL), with applications to liquid 3 He. Remarkably, the simple 
estimate (2) qualitatively applies in the degenerate case, since 
at T <C ep, where ep is the Fermi energy, the FL is a dilute 
gas of quasiparticles [4]. Moreover, Eq. (2) yields a singular 
temperature dependence 77 oc 1/T 2 . The latter originates en- 
tirely from that of the quasiparticle inelastic mean free path 
£ = UFTqp(T), with h/rqp oc T 2 /ep the scattering rate at the 
Fermi surface in 3d [4, 5], while v = vf is T-independent. 
Physically, the T — > divergence signifies the increasing re- 



sistance to the shear flow u(r) that causes distortions of the 
Fermi surface. 

The characteristic dependence 77 oc 1/T 2 , first estimated 
by Pomeranchuk in 1950 [1, 6] on dimensional grounds, was 
later confirmed by the calculations based on the FL kinetic 
equation [2, 3, 7]. It was also shown that the second viscosity 
£ is practically irrelevant, £ ~ (T/ep) 2r q [3]. Subsequent 3 He 
measurements [8] confirmed the relation 77T 2 — const. 

In the present work we consider the shear viscosity of 
the two-dimensional Fermi liquid. From a practical stand- 
point, the problem is relevant to a variety of Fermi sys- 
tems, ranging from the 2d electrons in heterostructures, to 
the trapped Fermi gases. Theoretically, the problem is com- 
pelling since in a 2d FL one generally expects nonanalytic 
energy- and temperature-dependence of response functions 
due to strong restrictions on the quasiparticle scattering [9]. 
Indeed, a naive estimate of the kind (2) should give a loga- 
rithmic suppression 77 ~ r qp ~ 1/T 2 ln(e_p/r) due to the 
well-known enhancement of the quasiparticle scattering rate 
Tqp 1 oc (T 2 /e F )ln(e F /T) [10-19]. Quite unexpectedly, we 
find that the 2d viscosity is enhancedby the square of the large 
logarithm ln(e p /T), as compared to that of the 3d FL: 
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The behavior (3) originates from the Cooper-channel pro- 
cesses in which the sum s = pi + P2 of the colliding mo- 
menta is much smaller than the Fermi momentum, s -C Vf- 
The logarithm arises from the corresponding singularity of the 
vertex part T oc F„/ln(pp/s) [4, 20], while the quasiparticle 
collisions at other angles appear to be less relevant. 

The temperature dependence (3) can manifest itself in the 
form of the interaction correction to the transport properties 
of 2d electron systems in the presence of a smooth disorder 
potential. The temperature-dependent hydrodynamic contri- 
butions to conductivity, a oc 1 /r)(T), caused by the resistance 
to the laminar flow of an electron liquid, have been anticipated 
since 1960's [21]. Recently the role of hydrodynamic modes 
in electron transport was revealed by observing the switch- 
ing from Knudsen to Poiseuille flow of the electrons in 2d 
wires [22, 23]. Assuming the laminar flow in the bulk 2d 
samples, the result (3) suggests that in d = 2 the hydrody- 
namic modes may lead to a singular temperature dependence 
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of the resistivity R xx oc In 2 T/T 2 . This prediction can be 
relevant for transport measurements in clean 2d heterostruc- 
tures in the metal-insulator transition regime: The interaction 
between the delocalized carriers may cause the apparently "in- 
sulating" correction to transport [24]. Viscosity may also limit 
the functionality of the FET devices suggested as a means of 
plasma waves generation [25]. Finally, hydrodynamic modes 
can play a role in the dynamics of the trapped atomic Fermi 
gases [26, 27]. The result (3) points at a possibility of mea- 
suring the Cooper-channel amplitude (defined below) as a 
functional of the repulsive interaction between the fermions. 

In what follows, we consider the 2d FL Boltzmann equa- 
tion, and draw a distinction between the scattering contribu- 
tions to the quasiparticle lifetime and to the viscosity. 

The Boltzmann equation. — The kinetic equation for the FL 



d t n + d r nd v e 
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St{n}, 



(4) 



where the quasiparticle energy e p 



,(o) 



Se, Se 



J f pp ' Sn p 'dT p ', is the functional of the distribution function 



n = rto(e(°)) + 5n, and f is the Landau function [dr p = 
2d 2 p/(2Tr~h) 2 ]. We assume the presence of a small nonuni- 
form velocity flow u(r) such that the equilibrium distribution 
function is 



no = f (e p 0) -pu-e F )/T 



f(x) = (e* + 1)- 



The linearized LHS of Eq. (4) is obtained by setting n : 
and keeping the flow gradients diUj in the second term: 



(5) 
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The first and the third terms of Eq. (4) almost cancel each 
other. Their difference ~ (T/e F ) 2 yields C ~ {T/e F ) 2 7j [3]. 

Quasiparticle collisions conserve energy and momentum. 
The collision term St {n} vanishes for the equilibrium dis- 
tribution no(e) = n — Sh, where e is the true quasiparticle 
energy, and the deviation Sh = Sn — (dno / de)Se includes the 
FL corrections. As usual, we identify the smooth part of Sh as 
Sh = —(dno/de)ip, where dno/de = — n (l — n n )/T. The 
resulting linearized collision term written in terms of ib has 
the form similar to that of a weakly-interacting Fermi gas: 

St{n} = f dT 2 dTyw(6)n 01 n 02 (l - n y){l - n 2') 

x6(e 1 + e 2 - e v - e 2 >)(ipi + ip 2 ~ ipv ~ tp2>) ■ (7) 

Here Pi + P2 = Pi' + P2' is assumed (see Fig. 1), 
9 = Z(pi,p2), and the scattering probability w(9) ~ 
(27r/7l)|2: 2 r(pi, p 2 — » pi',p2')| 2 is determined by the quasi- 
particle interaction vertex V and the quasiparticle weight z [4]. 

We look for a solution ip of the linearized Eq. (4) in the 
form of (6), V'(p) = ~ q(p)Y{P, u )^ where the angular part 



Y(p, u) = Pj 
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has the form of the second angular harmonic Y oc cos 2(<j> p — 
4>u), 4>p = &T g{Px + ip y ), with respect to the angle <p u = 

\ tan_1 dlZ-d 1 ^ defined b y the flow u ( x > V)- Substitut- 
ing ip into the collision term (7), we aim at the equation 
for the scalar deviation q that has a meaning of the vis- 
cous relaxation time. The integration measure drydTz = 
v F d£, 2 d^v{d6 /2-K)(d6y /27r), where £ = e — e F , and v F = 
m*/(irTi 2 ) is the 2d density of states. The energy delta- 
function sets the value for the angle 6\> — Z(Pi', Pi) via 

<5(ei +e 2 - e v - e 2/ ) = \de 2 . ld9 v \- 1 S(9 v - 9f) (9) 
with the Jacobian [10] 

dey/dOy = -2(eie 2 sin 2 6» + J 4) 1/2 ~ -2e F |sin0'| 
where 9' = Z(p 1 /,p 2 ') (Fig. 1), and 

A = - 66 = (6 - 60(6' - 6) • (ID 



(10) 



= In the last part of Eq. (10) we used cos 9' ~ (1 



A = A/e F 2 , that follows from P1P2 = Pi'P2'- The relation 
between 9 and 9' is illustrated in Fig. 2. 

The symmetry of the integrand in Eq. (7) with respect to 
9 — > —9 allows us to substitute Y 2 — > Yicos20, Yy — > 
Y\ cos 29y, and Y 2 i — > Y\ cos 2^2', and to cancel Y\ = 
y(pi, u) from the both sides of the Boltzmann equation. In- 
troducing the dimensionless energy variables x = £/T, as a 
result we arrive at the integral equation for q(x): 



y\T 2 



j(x 2 )f{-xv)f(-X2>)w{O)Qd0dx 2 dx; 
(2ir) 2 x 2e F |sin0'| 



(12) 

Here x 2 > = x\ + x 2 — xy, the integration in 9 is between 
and 7r since the particles are indistinguishable, and 

Q = q(xi) + q(xy)(cos26 — cos20y — cos26 2 ') 
= qi - qy + 2qy (sin 2 6y + sin 2 6 2 > - sin 2 9) . (13) 

The stress tensor = Jdr Shpide/dpj, combined with the 
definition (1), yields the viscosity in terms of q{x): 



q(x)dx 
(2 cosh f) 2 ' 



(14) 




FIG. 1: Quasiparticle collisions on a plane. At finite temperature 
T <C £f, the incoming pi, p2, and outgoing py, p 2 ' momenta 
can slightly differ in magnitude from p F . Left: a generic collision. 
Right: the Cooper-channel process, 9, 6' -C 1 and Oy ~ 6 2 > ~ 1. 
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FIG. 2: Left: Typical angular dependence of the quasiparticle scatter- 
ing probability w{6) that includes the In 2 \ir — 9\ Cooper-channel sin- 
gularity. Right: The relation between angles 9 and 9' illustrated, de- 
pending on sign (A), according to cos 9' = (1 — 4-) cos 9, 1^1 *^ 1. 



where N = vpep. 



and 



Above we used (piPj) 
(PiPjPkPi) = %(5ijSki + SikSji + SuSjk). 

Quasiparticle lifetime. — The solution of the integral equa- 
tion (12) is related to the problem of the 2d quasiparticle life- 
time T qp that is determined via substituting Q(x\, xy , x 2 , 9) 
by the isotropic term r qp (a;i). Below we will focus on this 
simpler problem first, to underscore the different roles played 
by the Cooper channel scattering processes in the inelastic 
lifetime and in the viscosity. 

To calculate r qp , one specifies the scattering probability 
w{9) entering Eq. (12). In Refs. [10-16, 18, 19] the scat- 
tering was assumed to be ^-independent, w — const. In this 
case the forward scattering (9 — 0) and the Cooper-channel 
(9 = 7r) contribute equally to the logarithmic singularity in 
^ via j;/ 2 ^ = J;"/ ^ * m(4|l|-V 2 ), \A\ « 1. 

Here 9 — for A > and sin 9 = V^A for A < 0. 

We now argue that the contributions of the forward scatter- 
ing and of the backscattering are in fact parametric ally differ- 
ent if one takes into account the Cooper ladder of diagrams 
that are logarithmically divergent to all orders in the inter- 
action. Summing the ladder causes a singularity in w(9 » 
7r;0y) ~ f |z 2 r( Pl , Pl ,;s)| 2 , 9y = Z( Pl , py), due to the 
renormalization of the Cooper-channel interaction vertex 



VpZ 



r(pi,pi-;s) 



E 

0,2,... 



cos m9y + cos m6y 

1,3,... 

(15) 

(Ref. [4], Sec. 20): In the leading-logarithm approximation, 
F m = f£ ] /[1 + Vm F^ In Ef], 1/ Vm = 2- 6 m0 , and F$ 
are the bare couplings. [Recent RG treatment shows [20] that, 
in general, even the asymptotic behavior F m = F m j ln(l/|0|) 
is characterized by non-universal parameters F rn ~ 1.] The 
above singularity suppresses the scattering probability w(8 ~ 
7r) ~ vo-xl In 2 where 8 = n - 9 < 1 (Fig. 2, left panel), 
and w n is dominated by the leading harmonic in Eq. (15). 
Thus the lifetime r qp is determined by the forward scattering, 
-0o MB)de _ w ^ ln ( 4 j^|— i/a), where we neglected the 

it/2 (F ( { ° )s ) 2 d\ni 



backscattering contribution ~ 
plitude w(0) is expressed in terms of the Landau parameters 



(l-F^l, 



The am- 



M = FP + F°oyj 2 , FP>°{9) = J2™=o Kf cosmO, as 



2tt 



4^(0) ~ — F 2 , 



F z 



r 



1 



1 + Fn° 



(16) 

Here we approximated the probability w(0) by its average 
over the Fermi surface (neglecting the 9 = it suppression) and 
averaged over the spin polarizations [28]. As a result [29], 



87T6F 



(17) 



The rate (17) is twice smaller than that of Refs. [15, 16, 19] as 
a consequence of the Cooper-channel renormalization [30]. 

The role of the Cooper channel in 2d. — The above calcu- 
lation demonstrates that out of the two singularities (at 9 = 
and 7r) of the integrand in Eq. (12), the former determines the 
lifetime r qp , while the latter is suppressed by the Cooper log- 
arithm. Below we argue that in the case of the viscosity, it 
is the 9 = 7r scattering that dominates, while other collision 
angles 9, including 9 = 0, are less relevant in T ftp. Indeed, 
for a generic 9, the energy and momentum conservation select 
the forward scattering, 9y ~ and 9 2 i ~ 9 (see Fig. 1). Due 
to the angular structure of Eq. (13), these contributions to the 
integrand are subdominant in the powers of T/ep since they 
originate from the deviations of the colliding momenta from 
Pf. This feature is specific to the 2d scattering. In contrast, in 
3d, the collision term for the viscosity is acquired from all the 
scattering angles 9 and cj>, where <fi is the angle between the 
planes defined by the incoming and outgoing momenta. (The 
3d collision term of Refs. [2, 3] indeed becomes small in T/ep 
for the in-plane scattering <j> = or it.) For the head-to-head 
collisions (9 = it), the phase volume for scattering rapidly 
increases (similar to the situation in BCS superconductivity), 
and there is no T/ep suppression of the integrand. 

Viscosity calculation. — The main difference of the colli- 
sion term (12) from that for the inelastic lifetime is the com- 
plicated angular structure in Eq. (13), which we will focus 
on below. In accord with the above discussion, the main 
contribution comes from the region where sin 2 9y ~ 1, in 
which case \9\, \9'\ <C 1, where 9' = tt — 9' (right panel of 
Fig. 1). In this limit sin 2 6y « sin 2 8 2 > > sin 2 9 « 9 2 , 
such that Q w q(x\) — q(xy) + 4q(xy) sin 2 9y. To esti- 
mate 9y in the limit \9\, \9'\ <C 1, we introduce the unit vec- 
tor s ~ (pi — P2,Pf8)/s in the direction of pi + P2 (here 
we chose the first coordinate axis parallel to pi), and the unit 
vector s' ~ (py ~ p 2 > ,Pf6')/ s that is rotated relative to s by 



o2 a>2 



9y sa fly. Heres 2 = {pi-p 2 )_+P F 8 2 _ = (py-P2>) 2 +P F 
Then cos6»i< ~ §§'. Since 9' 2 ~ 9 2 + A, it is clear that 
for sin 2 9y — 1, one needs either 9 2 < 9' 2 for A > 0, 
or 6' 2 < 

s/pf ~ max 



2 for A < 0. For such values of 9 and 9', 
9'\}, and we estimate 

2 



mm 



{{9/9'f, (9'/9) 2 } 



(18) 



[Dropping the (pi — p 2 )(py — p 2 >) terms in si' is justified 
since their contribution will be smaller in powers of T/ep af- 
ter the energy integration.] To perform the angular integration 



4 



in Eq. (12) we keep the lowest singlet and triplet harmon- 
ics in the series (15), and assume Fq ~ Fq/1d.(1/\6\) and 
F{ ~ Fl/\n{l/\e\) [20]. Then their contributions are 



L 



(F o s ) 2 d0sin 2 6 1 



! 2 In 2 VI 0' 



e=o 



o In 2 |0| \JA + sin' 

(Ff) 2 cos 2 9 v d9sin 2 6»i 

Jo 



In 2 VI' 

(^) 2 



'0 In |6»|V-4 + sin 2 6> 3 In 2 VI 



(19) 



(20) 



(here A > 0). For A < 0, the calculation is similar by in- 
terchanging 9 <-> 6', 6d6 = & d& ', yielding the above re- 
sults with A — ► \A\ under the logarithms. Thus the net an- 
gular contribution is jF 2 /\n 2 \Z\A\, where we introduced 
the spin-averaged dimensionless coupling \F 2 , with F 2 = 
(F S ) 2 + [F[f. We plug Eqs. (19) and (20) into Eq. (12), use 
the standard integrals Jdx-zdxy f{xi)f{— Xy)f(—X2') = 



f(-xi 



cosh(a:i/2) 
sh(x 1 ,/2) 



IC(xi — xy), where K(x) = xj (2sinh |), 



and fdxvdx 2 dx2> f(x2).f(~x v )f(-X2>) = y(-xi)(n 2 + 
x\), and obtain the integral equation 



1 



cosh^- 



2 +x 2 



\)y{xi)-j 



dxyiC{xi-xv)y{xv) (21) 



for the function y(xi) 
JC(x) =JC(x) x 2 



qQi) 



2r cosh(a;i/2) Ty/l+xj 
2 



1 



t 1 



In - j= = , where 

, (22) 



e F 



TVTT^ 



and ?i/r = T 2 /(47re F ) [29]. 

To solve Eq. (21), we note that the viscosity (14) is deter- 
mined by £ < T, which allows us to set x = under the loga- 
rithms, In 1/ VPT ^ ln(e F /TVl + x 2 ) w ln(e F /T), reduc- 
ing it to the standard problem [2, 7], whose solution yields 



T„ = 



TO 



21n^ 



fl 4a v-^ 



2n+ 1 



1 



n 2 (n + l) 2 n(n + 1) — a 



(23) 

where a = 2 [l - (F^/Fo) 2 / ln 3 (e F /T)] . Since 2 - a < 1, 
from Eq. (14) we obtain the result (3). As advertized, the vis- 
cosity is fully determined by the quasiparticle scattering prob- 
ability v 2 F w ~ x|^ 2 / ln 2 (e F /T) in the Cooper channel. 

In summary, we found a nonanalytic temperature behavior 
of the 2d Fermi liquid viscosity and related it to the logarith- 
mic singularity of the quasiparticle scattering amplitude in the 
Cooper channel. The ratio between the viscous and the inelas- 
tic scattering times is enhanced by the factor <~ ln 3 (e F /T) 
due to the restrictions on the 2d quasiparticle scattering. 
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